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Summary 

The initial design and demonstration of an Intelligent Flight Propulsion and Control System 
(IFPCS) is documented. The design is based on the implementation of a nonlinear adap- 
tive flight control architecture. This initial design of the IFPCS enhances flight safety by 
using propulsion sources to provide redundancy in flight control. The IFPCS enhances the 
conventional gain scheduled approach in significant ways; 

• The IFPCS provides a back up flight control system that results in consistent responses 
over a wide range of unanticipated failures. 

• The IFPCS is applicable to a variety of aircraft models without redesign and, 

• significantly reduces the laborious research & design necessary in a gain scheduled 
approach. 

The control augmentation is detailed within an approximate Input-Output Linearization 
setting. The availability of propulsion only provides two control inputs, symmetric and 
differential thrust. 

Earlier Propulsion Control Augmentation (PCA) work performed by NASA provided 
for a trajectory controller with pilot command input of glidepath and heading. This work 
is aimed at demonstrating the flexibility of the IFPCS in providing consistency in flying 
qualities under a variety of failure scenarios. This report documents the initial design phase 
where propulsion only is used. Results confirm that the engine dynamics and associated 
hard nonlinearities result in poor handling qualities at best. However, as demonstrated in 
simulation, the IFPCS is capable of results similar to the gain scheduled designs of the NASA 
PCA work. The IFPCS design uses crude estimates of aircraft behaviour. 

The adaptive control architecture demonstrates robust stability and provides robust per- 
formance. In this work, robust stability means that all states, errors, and adaptive parame- 
ters remain bounded under a wide class of uncertainties and input and output disturbances. 
Robust performance is measured in the quality of the tracking. The results demonstrate the 
flexibility of the IFPCS architecture and the ability to provide robust performance under a 
broad range of uncertainty. Robust stability is proved using Lyapunov like analysis. 

Future development of the IFPCS will include integration of conventional control surfaces 
with the use of propulsion augmentation, and utilization of available lift and drag devices, to 
demonstrate adaptive control capability under a greater variety of failure scenarios. Further 
work will specifically address the effects of actuator saturation. 
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1 Introduction 

The objective of this work is the design and demonstration of an intelligent model inversion 
controller architecture that provides control augmentation using propulsive forces only. The 
controller does not require extensive tuning and is applicable to a variety of aircraft without 
redesign. It provides backup flight control and maintains consistent response over a wide 
range of uncertainties and disturances, including unanticipated failures. 

This report describes the lateral and longitudinal applications. The applications aim to 
provide redundancy to the conventional control surfaces. The possibilities of using propulsive 
control to provide reasonable handling are explored. The application of propulsive control 
is limited by throttle saturations in both rate and position. Furthermore, the bandwidth of 
the augmentation is governed by the engine dynamics. 

One of the objectives is to implement the controller architecture on the “Stone-Soup- 
Simulator” (SSS) as used by NASA Ames. This provides the proper environment for in- 
tegration onto the Advanced Flight Control Simulator (AFCS) used at Ames. The code 
associated with the SSS in its current version does not provide a convenient environment for 
development and analyses. Therefore, a similar nonlinear dynamic model representative of a 
civilian transport aircraft was used for the analyses in this report. The results of integration 
into the SSS code are included as a seperate chapter. 

1.1 The NASA PCA program 

Following the DC-10 accident at Sioux City, Iowa in 1989, the NTSB recommended; 

Encourage research and development of backup flight control systems for newly 
certified wide-body airplanes that utilize an alternate source of motive power 
separate from that source used for the for the conventional control system [1]. 

In general there is no satisfactory method onboard the aircraft for effectively controlling 
the flight path with a disabled primary flight control system. Manual throttle control of 
the engines is extremely difficult due to the dynamic response of the aircraft in this situ- 
ation. NASA Dryden Flight Research Center (DFRC) investigated the use of propulsion 
for emergency flight control in a variety of aircraft. These investigations demonstrated that 
throttle control of engines alone can be used to augment or replace the aircraft primary' 
control system to safely land the aircraft. DFRC used specific control laws for bank angle 
and flightpath angle in the flight control computer, driving the engines in response to pilot 
input commands. An example as applied to a Boeing 747-400 is presented in Fig. 1.1, taken 
from Ref. [2]. The gains in Fig. 1.1 are scheduled with respect to 

• velocity, 

• altitude, 

• c.g. location, 

• operating mode: 


1 



Figure 1.1: PCA conventional control law block diagram. Gains are scheduled with: ve- 
locity, altitude, c.g. location, operating mode, configuration, nature of failure, location of 
operational engines. 
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- unusual atitude recovery, 

- cruise flight, 

- ILS-approach, 

• configuration: 

- flaps, 

- landing gear, 

• nature of failure: 

- jammed vs. 

- floating control surfaces 


and, 

• number and location of engines available for control. 

The experience of the designers, as well as access to accurate simulation models, allowed 
for good gain scheduled designs based on classic linear control techniques. However, gain 
scheduling with respect to the nature of failures significantly limits the safety that may be 
provided by the controller structure. It requires not only that the operator knows which 
failure condition is current, but also anticipation of the particular failure. 

The applicability of a back-up flight control system will be greatly 
increased if it can provide desirable responses over a wide range of 
unanticipated failures. 

In flight failures may include partial loss offlifting surfaces, or asymmetricaly jammed control 
surfaces. Furthermore, the design benefits from avoiding the laborious research necessary for 
a gain scheduling design approach. It can relieve the reliance on accurate simulation models, 
which are not always available. 


1.2 Nonlinear Adaptive Propulsion Control 

The work presented in this report aims at the implementation of a nonlinear adaptive flight 
control architecture that provides a flight safety enhancement by using the propulsive sources 
as redundant control. The nonlinear adaptive architecture improves upon the conventional 
gain scheduled approach in design of a PCA, by allowing for a large variety of in flight failures 
with a minimum of gain scheduling. 

A rigorous theoretical development is provided. The theory combines recent develop- 
ments in approximate linearization and nonlinear adaptive control. The development is 
based on a proof using Lyapunov like stability analysis. The proof is not included in this 
report but can be found in Ref [3], 
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1.3 Transition of Results 


For syntheses, we used simulation involving a linear model of the MD-11, as well as a 
nonlinear model of a generic airliner in Mat Lab-environment. The linear representation of the 
MD-11 is from Ref. [4]. The nonlinear aircraft model is based on the European Research Civil 
Aircraft Model (RCAM) [5]. This aircraft model represents an initiative designed to promote 
the use of robust flight control design techniques in the European aeronautical industry. The 
RCAM includes actuator and engine nonlinearities, as well as nonlinear aerodynamics. The 
aircraft weights 220,000-331,000 lbs. It has two engines underneath its main wings, together 
capable of providing a thrust to weight ratio of approximately 0.35. The engines have no 
inclination angle, and are about 6 ft below the eg, and 26 ft outside the plane of symmetry. 
The aircraft is configured for landing. With these characteristics, the RCAM offered a 
convenient basis from which to investigate various aspects of the PCA application. 

The results obtained within the Matlab environment were subsequently adapted and im- 
plemented within the portable version of the miniACFS-code of NASA Ames. This portable 
version is referred to as Stone Soup Simulator (SSS). Due to the limited access to information 
pertaining to the SSS, it was not practical to obtain reliable stability & control derivatives 
of the aircraft represented in the code. The nature of the controller presented in this report 
allows for good performance, under mild assumptions, and using a crude inversion model. 
However, prudent design practice suggests the use of the most accurate data available. Ini- 
tially, a Jacobian linearization of the RCAM was used for model inversion control design for 
this code. A better approximation was later obtained by utilizing Boeing-747 stability & 
control derivatives [6] within the linear model inversion. 

It is expected that the controller architecture will be further evaluated in piloted simu- 
lated flight at NASA Ames research facilities in the second phase of this project. 
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2 Input Output Linearization with 
Neural Network Augmentation 

2.1 Introduction 

The controller architecture proposed for this part of the research is based on Feedback Lin- 
earization (FBL). FBL has attracted a lot of attention in the last decades. Fundamental 
contributions have been provided in several text-books, including [7, 8, 9]. This method is 
fundamentally different from conventional linearization, as it is an exact state transformation 
and uses state feedback, rather than local linearization of the dynamics. FBL algebraically 
transforms the nonlinear dynamic system into a linear one. Linear control techniques can 
then be applied to the transformed system. 

In exact FBL, the outputs are differentiated with respect to time (or the independent vari- 
able) until the first explicit dependence on the control. The number of required derivatives, 
r, for a given output variable is its relative degree. However, sometimes the first appearance 
of the control is weak or negligable for reasonable values of the control variable. The FBL 
theory can still be used in the form of Approximate Feedback Linearization, formalized by 
Hauser [10], and extended to non-control affine systems in Ref [11]. 

Approximate FBL combines the advantages of a conventional linearized approximate 
model of the aircraft dynamics, within the FBL paradigm of the full nonlinear system. 
Using an approximate model, rather than exact algebraic expressions results in an inversion 
error. We augment the approximate FBL with an adaptive ’learning-while-controlling’ neural 
network (NN) to compensate for this error. 


2.2 Approximate Feedback Linearization and Bounded 
Tracking 

Consider the aircraft dynamics represented in the general form 

x = f{x, u) (2.1) 

y = h(x, u) (2.2) 

We wish to design a tracking controller for a smooth and bounded trajectory with bounded 
derivatives of sufficient order. Suppose that in the i 01 element of the output the components 
of are all small. We reorganize those elements of the output function as 


y» = 4>t, o(x) + expifiix, u) (2.3) 

so that erpifi « <f>i t o< With this representation, the output may be differentiated successively 
until the control effect is significant. Each output y* is differentiated r { times. Starting with 
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Eqn. (2.3), 


i/i 

= 4>i, i(x) 

+ 

eA, i(x,u,u) 

• 

ss : 

+ 

* 


= 4>i, r,-l(x) 

+ 

u, • 


= ^(x.u) 

+ 

erp i>Ti {x, u, u,-*0 


(2.4) 


We say that the effect of the control is significant when = 0(<f>). At each step, the terms 
that are significant are included in <f> and those that are not are included in eip. This means 
that V'ij includes, 


• the time derivative of tpij-i, 


• small terms from the time derivative of 


Definition 2.1 A system of the form 

x = / (x, u) x € 1R D , u, y € !R m , 2 

y = $ 0 (x) + e^o(y,u) 1 ' 

is said to have a well-defined approximate local vector relative degree [ri,r 2 , • • • , r m ] T in the 
neighbourhood (fi* x Q. u ) of a point (x 0 , u 0 ) if the Jacobian is nonsingular Vx € fi* and 
Vu € fi u , where 

*(x,u) = [*,,„••• (2.6) 

and, where r { is the least order time derivative ofyi on which the effect of control is significant. 


Theorem 2.1 Approximate input-ouput linearizability 

For a system of the form (2.5) that has a well-defined approximate relative degree [r l5 r 2 , • • • , r m ] T , 
there exists a state- dependent control transformation to a pseudo-control 

!/• = [■,; „• ... C] r = <f ( x ,u) (2.7) 


that transforms the system into m decoupled approximately linear systems 

yS r,) = + e^i >ri (x, u, u, • • ■) 

Pm m) = ^ + U, U, *) 


( 2 . 8 ) 


A proof is presented in [11]. 

Consider a system of the form (2.5) that is approximately linearizable. The total approx- 
imate relative degree is defined as 


r 



«=i 


(2.9) 
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If f = n, then the closed loop can be interpreted as approximately linearized from input- 
to-state by transformation to state- variables £, where 

(i = [<Mx) ••• *,,-.(>)!' (210) 

( = [ff£-£] T (2ii) 

If f < n, then the nonlinear system (2.5) can be transformed into a so-called normal form. 
The normal form of the system may be represented as 

ti = [6,2 Zi, Ti $(£,^,u) + e$(£,/i,u,iv-0f (2.12) 

ii = ©(*,/!) (2.13) 

With f < n, there are internal dynamics, which must be examined for stability. In general, 
these internal dynamics may depend nonlinearly on the normal outputs, £. Therefore, we 
consider the so-called zero-dynamics , which are the internal dynamics when the output y, 
and all its time derivatives are maintained at zero by the control u, which implies £ = 0. 
Thus the control u is such that the zero-dynamics are formally defined by 

U = 0(0, /i) (2.14) 

Remark 2.1 It is not necessary to obtain the explicit expression for fi. We may also analyse 
the internal dynamics by completing the state transformation Eqn (2.12), with n-f vector 
fields, and using the expression for u(x) that provides for £ = 0. 

Linearized Model Inversion Design 

For convenience define £ r , 


tr = [2/i n) ••• y£ m) ] T (2.15) 

The approximate linearizing transformation of definition 2.1 and theorem 2.1 may be sum- 
marized as 


where, similarly to Eqn 2.6, 


ir 

= v* + e^(x,u) 

(2.16) 

1/* 

= #(x,u) 

(2.17) 

*(x,u) 

= V’m.rrn ] 

(218) 


If$(x,u) is invertible with respect to the control, we may express the original control input 
in terms of the pseudo control 


U = $ 1 (x, v) (2.19) 

Considering the nonlinear nature that we wish to allow for the plant, the inverse of $(x, u) 
may be particularly complex are impractical to obtain. If the condition of Theorem 2.1 is 
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satisfied, then the Jacobian $ 0u = js£| is nonsingular for a set (xo,uo) € (Q* x Q u ), 

m lxo,tio 

and a linearized form of 2.19 may be constructed as 

u = uo + $ou {v ~ $o - $ox ( x - x o)} ( 2 - 20 ) 

where uo and x 0 are trim values of the plant, $o = $( x o>Uo)> and $ox is the Jacobian 

with respect to the state at (xo, Uo)- Using Eqn 2.20 for the control law implementation will 

result in an inversion error. If we define this inversion error as 

A(xo,Uo,x,i/) = $(x,u) - {$o + ^o*( x_x o) + ^oti(u-Uo)} (2.21) 

Remark 2.2 Note that Eqn (2.20) may in fact he implemented with approximations for the 
trim values, the Jacobians, and the nominal operating point 

u = Uo + l>ou {i/ - $o ~ So* ( x “ *o)} ( 2 - 22 ) 

The approximations should be designed ‘reasonably’ accurate to avoid unnecessary control 
effort. Specifically, the sign of the elements of the Jacobian should be known, 

spn(^|y) = sgni^lij) (2.23) 

Using Eqn (2.21), Eqn (2.17) can be expressed as 

v = $(x,u) - A(x 0 ,u 0 , x ,i') (2.24) 

and the closed loop as 

£ r - v + A £ (x 0 ,uo, x ,v) (2.25) 

where 

A e (x 0 ,uo, x ,*') = A(x 0 ,Uo,x,i/) +e^(x,u) (2.26) 


2.3 Online Neural Network Augmentation 

In the first phase of this work we use a separate neural network (NN) for each channel. Thus, 
the NNs have a single output. The outline provided here is stated for NNs with multiple 
outputs, allowing a single NN for multiple channels. In future work we will consider the 
extension of the application to include a single NN for multiple channels [3]. 

We design the pseudo control signal as 

V = Vpd - v ad - v T (2.27) 

where Vad is an adaptive contribution from an online approximator to be detailed next, and 
v T is a robustifying signal. Construct i ' pd = [ v p d x • • • Vpdr* ] T j where if i — 1, • • • ,m 

(2.28) 

P= 1 
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where the coefficients k ip chosen such that the error dynamics y\' = v v ^ are exponentially 
stable. Define A { G R(n-»>*(*»-i> and b< G R.(n- l )xct) ^ 


f 0 

l 

0 

0 ••• 

0 ' 


' 0 ' 

0 

0 

1 

0 ••• 

0 

bj = 

0 


ki 2 


• • . 



. 1 . 


Ai = 


Define = y[ p ^ — y^ p \ and also define 

rp 

ei = [n y (r ‘ _l) ] 

With these definitions, the elements of Eqn (2.25), which are 

y (ri) *= Vi + A ti 

may be represented in state space canonical form as 

dj = A&i + bj + v r — Ae,) 


(2.29) 


(2.30) 


(2.31) 


(2.32) 


2.4 Linear Neural Network Structure 

This section derives a network adaptation law that guarantees boundedness of the tracking 
error and of the network weights. Consider a linear-in-the-parameters NN approximation of 
a possibly nonlinear inversion error 


W*'p = f-e (2.33) 

where 0 < j|e|| < E. The vector e is referred to as the NN reconstruction error, or residual 

error. W* is an ni x ^-matrix of constant ’ideal’ parameter values that, in some sense, 

minimize ||e||. These values are not necessarily unique. (3 is an ni-vector consisting of basis 
functions. Here n? gives the number of NN outputs. 

For a linear-in-the-parameters network 

u AD = W'P(x,u) (2.34) 

where W is an n\ x ^-estimation of W*, and x is made up of selected elements of the state 
vector. The selection of the elements of x is done through careful assessment of the inversion 
error [12]. Since the input to the NN includes the pseudo control signal, v , which in turn 
is a function of the output of the NN, u A o , a fixed point problem occurs, v is input to the 
NN through a squashing function. This insures that at least one fixed point solution exists. 
The error dynamics can be written with ri 2 = 1 as 

e = Ae + b(v AD - f) (2.35) 
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with 


A = 


0 1 
-K P -K d 


and e = [ / q q ]' and b = [0 1 ]'. Define the weight error as W = W — W*. Rewrite 
as _____ 

e = Ae + bW'0 + b (W'0 - f) 

The update law of the NN weights is designed as 


w = w = -rc/3 

T = V > 0 and where C is a scalar filtered error terra, 

C = e'Pb 

with 

P = 


r . _i_ _j_ 

K P T 1K d 2Kp 

1 14Kf 


2K P 


2K p K d 



'Pi 1 

P ) 12 1 


P\2 



(2.36) 
(2.35) 

(2.37) 

(2.38) 

(2.39) 

(2.40) 


Remark 2.3 The implementation as presented provides for ACAH in the longitudinal chan- 
nel, i.e. a design with second order tracking error dynamics. This design may be generalized 
where 

e = [x (n) i (n “ 1) (2.41) 

and b = [ 0 • ■ • 0 1 ]' a size-n vector, and P an n x n-matrix that solves 

A' n P + PAn — 2 I n = 0, with An in canonical form similar to (2.36) but expanded to order n. 


Remark 2.4 Notice that the closed loop system defined by Eqns (2.37), (2.38), and (2.39) 
may also be generalized to n 2 > 1. This would be the case if one NN is used in a controller 
architecture that integrates the three control channels. In that case we can, for instance, 
organize 

e = [x'x / ]' (2-42) 

where x and x are size ^-vectors. A will then be designed a 2n 2 x 2 n 2 block matrix, 


A = 


0 ^4 12 

A21 A22 


(2.43) 


with A 12 — I ni and 


A 21 = — diag{Kp 1 ,Kp 2> - •• ,K Pni }, and 
A 22 = - diag{K Di , K Di , • ■ • , K Dni } 

where the index corresponds to the elements of x. Furthermore, b= [ 0 I n2 ]' and 

W = W = -T0E' (2.44) 


where E=[ Ci C2 — Crrt]'- 
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Neural Network Input Design 

The NN can consist of any linearly parameterized feedforward structure that is capable of ap- 
proximately reconstructing the inversion error. Ref [13] uses Radial Basis Functions (RBFs) 
because these functions are universal approximators even when the network is linearly pa- 
rameterized. However, it is well know that RBFs are very poor at interpolation between 
their design centers, and a large number of such basis functions are needed for networks 
with multi- dimensional input vectors. In Ref [14], RBFs were used to capture variations in 
Mach number, because in the trans-sonic region, these variations are difficult to represent 
by polynomial functions. In the current implementation, a single-layer sigma-pi network is 
used. The inputs to the network consist of the state variables, the pseudo control and a bias 
term. Fig 2.1 shows a general depiction of a sigma-pi network. The values represent the 
weights associated with a nested kronecker product of input signal categories, and therefore 
they are (binary) constants. The values are the variable network weights. 


kronecker 

product 



Figure 2.1: Neural network structure, designed for longitudinal channel. 

The input /output map of the NN for the longitudinal channel is represented as 

vad = W'(3(x, Vj, vp) (2.45) 

where Vg, and vp are respectively the pitch channel and roll channel elements of the pseudo- 
control v. The vector of basis functions, /9 is akin to the regressor vector in adaptive control 
texts. The basis functions are made up from a sufficiently rich set of functions so that the 
inversion error can be accurately reconstructed at the network output. The basis functions 
are constructed by grouping normalized inputs into three categories. The first category is 


11 





used to model the inversion error effects due to changes in airspeed, since the stability and 
control derivatives are strongly dependent on dynamic pressure 

C\ : {b u V T ,V>} (2.46) 

Where bi indicates a bias term. In allowing the plant to be nonlinear and uncertain in the 
control as well as in the states, the inversion error is a function of both state variables, Xi, 
and pseudo control, v. These, and a bias, are therefore contained in the second category 

C 2 : [bz, Xi, Xj, •••, v} (2.47) 

The third category is used to approximate higher order effects. These may e.g. be due to 
the transformation between the body frame and the inertial frame 

C 3 : {b 3 , Xi } (2.48) 


Finally, the vector of basis functions is composed of combinations of the elements of C\, C 2 , 
and C 3 by means of the kronecker product 


where, 


/3 = kron (kron(Ci, C 2 ), C 3 ) 

(2.49) 

kron(x, y) = [xij/iXij/ 2 Xm!/n]' 

(2.50) 
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3 Controller Design 

This section describes the approximately linearizing controller design for rate command 
systems in pitch and roll, as well as a yaw rate command system as alternative to roll rate 
command. It is assumed that the only available control consists of thrust from the outboard 
engines. 


3.1 Architecture 


[A P)c „ 


linear 
* >-► controller 


AD\ 


v p 


Pc 

command 

filters 4c 


Qc linear 

controller 





^Thl 

>-l 

J 7Vott/« 

" 

^Thl 


W 



— V AD2 


\p>p] 


Aircraft 


Figure 3.1: Architecture for angular rate commands using propulsive control. 


Figure 3.1 displays the architecture for the propulsive-control-only setup for rate-command- 
attitude-hold (RCAH) in pitch and roll. The inversion with respect to the throttles is based 
on the symmetric and asymmetric (or differential) application of thrust. These create re- 
spectively a moment about the body pitch and yaw axes. The roll rate is therefore indirectly 
controlled by the application of asymmetric thrust. 

Let the aircaft state be 

x = [pqruvw] T (3.1) 

We may represent the rigid body aircraft dynamics as 

x = [/i(x) / 2 (x) / 3 (x) / 4 (x) / 5 (x) / 6 (x) ] T (3.2) 


With the following inertia terms, 

h 


C3 = 


U, - il. 


J ^4 


hz _ 1 h 

TT _ J2 ’ c? / ,C9 T T — P 

1 x 1 z 1 xz J v 1 X 1 Z *xz 


(3.3) 
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let T, A, and, ft represent the aerodynamic moments, constructed as 


(3.4) 

(3.5) 

(3.6) 


r = c 3 l + C4N 
A = c 7 M 
ft = c 4 L + C 9 N 


where L, M, and N represent the usual moments about the body fixed axes. Then 


r/i 1 


(cir + c 2 p)q 

+ 

r ' 


0 

c 4 ' 

f a 


ChP r- -c#(p ? ■ -r 2 ) 

+ 

A 


Cl 

0 

h 


{c s p-C2r)q 

+ 

ft 

+ 

0 

c 9 

f 4 


Fx/m 



l/(z e m) 

0 

h 


F y /m 




0 

0 

. h . 


Fz/m 




0 

0 . 


M e 

. N e 


(3.7) 


where Me and N E represent respectively the moments due to symmetric and asymmet- 
ric thrust. We choose pitch and roll-rate as outputs for input-output linearization using 
symmetric and asymmetric thrust as input. 



(3.8) 


The effect of asymmetric thrust on roll rate is mostly a secondary effect. To obtain significant 
control of roll-rate using asymmetric thrust, we neglect the effect of the inertial coupling c 4 
in Eqn 3.7, and therefore consider the second derivative of y\. 


i/i 


dfi . dfi . df x . „ . 

a* p+ W ,+ ft rr + r ’“ 


The effect of thrust appears in the expression for y x as 


(3.9) 


Vi ~ r rC9 N e + (c^ + c 2 c 4 ) qN E + (cjr + c^p) M E (3.10) 

In the pitch channel, M e has a significant effect on the first derivative 

y 2 = c b pr - ce (p 2 - r 2 ) + A + c7 M E (3.11) 

For the following analyses we write the system in normal form [7]. Therefore define 
£ = [p P q] T a^d let 



Then, the input-output transmission may be presented as 




^(x) + ff(x) 


M b 

N e 


(3.12) 


(3.13) 
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Me 

N e 


Let represent an approximation, then 

£ r » £(x) + £(x) 

The feedback linearizing transformation is 

& = 

= ?M + S(x) 

Me and Ne may be obtained in terms of the pseudo control by inversion of 3.16. 


(3.14) 


L v i 


v i 


Me 

N e 


(3.15) 

(3.16) 


M e 

N e 


The inverse of Q is 


where det(^) 


Q- x = 


0 


x 

Cl 


C7 (cir+cap)c 7 

L — det(fil) -det(C) 


det(£) = - (( c x q + r r )cg + (c2? + r p )c 4 ) c 7 
Matrix Q becomes singular for an unlikely pitch rate, 


q = _£®ElZL^iEe > 0.5 rad/sec 

CjCg + C2C4 


(3.17) 

(3.18) 

(3.19) 

(3.20) 


In order to analyze the remaining dynamics when £ is controlled, we consider the decou- 
pled linearized modes in the lateral and longitudinal channels. If we choose the design based 
on these modes, Q reduces to 

0 1 1 (3.21) 


£>-‘ = 


1 CT 
1 0 


L r r c9+r p c4 

The symmetric and asymmetric throttle deflections, resp. s5th, and are obtained as 


s6th 

a5th 


= e--‘{[^]-^ (x) } 


where 


'th 


G' 1 = . 

2T r n ith+ 2T p r ith 

The closed loop stability analyses are contained in section 4. 


0 

1 


2A, 

0 


(3.22) 


(3.23) 
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3.2 Construction of the throttle positions 

From the signals a5 th and s5 th the throttle positions can be constructed as 

6*1 = |(4 + 4) (3-24) 

$th2 — g ( 5< ^ — °^) (3.25) 

3.3 Construction of the Pseudo Signals 

The transformation 3.15-3.17 allows the design of desired pseudo signals. The pseudo signal 
for Rate-Command-Attitude-Hold in the pitch channel can be constructed as 

Vq = q c + on q + a 0 qdt (3.26) 

where q = q c — q, and where q c , q c are obtained using a 1'* order command filter 3.1. In the 
roll channel, we require a 2 nd order command filter to generate p c , p c , and p c 

vp — pc J r Ct2p 4- Qip + a 0 £pdt (3.27) 


16 


4 Stability Analyses 

4.1 Introduction 

Most air-transport aircraft are designed with attitude and speed stability, but may display 
roll divergence. In the longitudinal direction, an aircraft is trimmed for a certain speed. 
Given a constant aircraft configuration, each combination of control surface positions will 
result in a particular speed. Speed stability implies that, given an external disturbance, the 
aircraft will eventually return to this unique speed. 

The trim point, as well as the stability characteristics, may differ between the cases of 
‘frozen’ or ‘jammed’ control surfaces and ‘free floating’ surfaces. In working with the SSS 
version of the miniACFS, there is no provision for simulating free floating control surfaces. 
Therefore, we used the trim positions provided by the SSS code and ’failure of the control 
surfaces’ simply implies that the control surface positions remain frozen from the time of the 
failure. 

This chapter provides the analyses of the PCA control laws and the associated zero dy- 
namics for the different channels. We apply input-output linearization techniques to provide 
rate command attitude hold (RCAH) control, and investigate its feasibility using propulsion 
only. 

The stability analyses are performed based on a Jacobian linearized representation of the 
aircraft in a landing configuration. This approach will provide for a straightforward means 
to analyze the zero dynamics. 


4.2 Roll channel Augmentation 

4.2.1 Linear analysis of zero dynamics 

The Jacobian linearized form of the lateral dynamics, without the engine dynamics, may be 
represented as 


' p ■ 


r 


V 


A. 

= A 


ft 


r r 

fl r 
r'p Y r 
1 tan do 
i + b a5th 


r„ 

flt) 

Yv Y* 
0 0 



' P ' 


2r ^ 1 


r 

4- 



V 

i 

2Y ith 

. 

A. 


0 


a5th (4.1) 

(4.2) 

The effect of asymmetric control input in roll is T Sth /Q Slh « 1. Following the definitions 
in Eqns 2.11, and 2.12, we obtain for the lateral channels 

e = [ppf (4.3) 

The zero dynamics may be analyzed by using £ = 0 and following Remark 2.1. By neglecting 
the effect of and Ys th we may approximately represent the zero dynamics by 
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[* j 



v _ LtH 

i v Tr 


Ys- 


r UV, 


— tan 6q — tan 8q J 



r*i 


i 


(4.4) 
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Example 4.1 For the MD-11 as presented in Ref [4], Eqn (4.4) has eigenvalues at —1.415 
and .005. The slow RHP pole is due to the roll divergence mode, and reflects the fact 
that, given a roll rate command system, the pilot provides the integration of this mode 
when tracking a desired bank angle. The effect of neglecting and Yf th , is in the 
order of magnitude of l.Oe — 3. 

4.2.2 Analysis of output tracking dynamics 

If we use a linear controller design for the pseudo control up 

vp = Pc + ct 2 (jp c - P) + ctiipc - p) + cto{<f> c - <f>) 

Then if we control p c = p c = p c = <f> c = 0 we obtain 

vp — - (a 2 p + ai p 4- Oq4> ) 

= - [o 2 r p + Q 2 r r a^r,, a 0 o ] x< at (4.5) 

where x/ 0( = [pru^f. We may check the closed loop modes by examining the eigenval- 
ues and corresponding eigenvectors. Choose the roll rate as the output of interest, 


C = [1 0 0 0 0] 

(4.6) 

Then CA represents the 1 st row of A, and 

p — CAx =*> 

(4.7) 

p = CAx 

(4.8) 

= CA { A x + b a5 t fc} 

(4.9) 

The linearizing transformation is represented by 

P = vp 

(4.10) 

up = CA {A x + b 0(5#,} 

(4.11) 


The first term A x will serve to cancel the natural modes of the aircraft. For various reasons, 
we may not want to use the complete linear model of the aircraft in the inversion. For 
example, in case of non-minimum phase zeros in the system, an inverting controller would 
be unstable. Also, we may choose not to cancel beneficial cross-coupling in the system. 
Therefore, allow for approximate linearization [15, 10], expressed as 

up = CA {Ax + baij/J (4.12) 

The inverse control transformation Eqn (4.12) is 

a5tk = (CAb) -1 {up - CAAxj 
= — kj x 
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(4.13) 

(4.14) 



In terms of estimates of aerodynamic derivatives, this becomes 


n e = ^ ^ - ( r „ r p + r r n p + r v Y p /m) P 

— (r p r r + r r n r + t v y t /tti) t 

— (r p r,; + T r Q v + T v Y v /m) v } 

= r iP r 2^ r 3^ } (4.15) 

1 r C9 

The lateral modes in the closed loop are then represented by 

A c = A - bki (4.16) 

We may analyze the closed loop dynamics about the equilibrium by checking if A c is Hurwitz. 

Example 4.2 For the MD-11, with £ = 0.8, =^0.5 and a 0 = 0, ax = u^, 

and ct 2 = 2(,oj n , applying ‘exact’ inversion by using A = A in Eqn 4.13, we obtain 
eigenvalues of Ac as 

-1.645 - 0.379 ± 0.273i 0.0019 (4.17) 

The unstable pole is due to the slightly non-minimum phase zero of the open loop 
aircraft. To avoid the singularity this causes in model inversion control architecture we 
can now use an approximate inversion by using A in Eqn (4.13), where A = >l| rr _ 0 . 
The resulting closed loop then gives poles located at 

-2.094 - 0.130 ± 0.209i -0.548 (4.18) 


4.2.3 Effect of engine dynamics 

The preceding results do not consider the engine dynamics. The engines are modeled as 
first order systems, with a time constant r eng . Let ay t h be the asymmetric output of the 
approximate engine model in terms of throttle position, 

dith = K en g o,Xth + K en g(i5th c (4.19) 

ayth = ax th (4.20) 

where K eng = T~„ g . Define x' = [x ax^ ] T - A linear representation of the closed loop then 
may be constructed as 


x = 


A b 
0---0 -K eng 


The closed loop is then represented by 


A' 


A b 

0---0 -A'ens 


x' + [ 0 • • • 0 K eng ] aS thc 

-[0---0 tf„„] 7 '[k 1 o] 


(4.21) 

(4.22) 
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The resulting eigenvectors & values show the effect of the engine dynamics on the feedback 
linearization. 


a: 
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eng 


r r c© 
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r r 
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0 0 2J 1 ,* 

0 0 2Q ith 
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T7c# Q ° ® ~Keng 


where 


A = A + aj -I- a 2 r p 
A = r 2 + a 2 r r 
A = r 3 + a 2 r„ 

with Ti i2i 3 defined in Eqn 4.15. 

Example 4.3 For the MD-11 [4], with £ = 0.8, w n = 0.5 and c*o = 0, c*i = wj(, and 
a 2 = 2£a; n , and using the approximate inversion A tat = ^ io tlr r =o Eqn (4.13), we 
obtain eigenvalues of A' c as 


-1.049 ±1.828* - 0.112 ±0.195i -0.616 (4.23) 


We may analyze the effect of the engine dynamics in a simple root locus representation, 
Fig 4.1. This will provide a straightforward tool to aid in choice of the gains by setting the 
ratio of K P to K D , thus providing a root locus by varying K D . 


Linear controller design 

* 


K is 


K p 


+i 


\s 




v iwtBian 
tewrolkw 




Arcnft 


Figure 4.1: Roll channel setup for design of the tracking error dynamics. 
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4.3 Longitudinal channel 

The longitudinal modes, without engine dynamics, are approximated by 


' 9 ' 


a 9 

A u Ay, 

0 * 
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to 

p* 

. 6 . 


1 
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0 J 


= Ax + b sSth 


(4.24) 

(4.25) 


where s5 t h represents symmetric throttle deflection, i.e. the sum effect of both throttle 
deflections, where &6 t hi = A5 t h2, and 2 Ag th sS t h = c-jMe with Me the moment about the 
body Y-axis due to engine thrust. 

Let Ai >: represent the first row of A, and bi the first element of b. The feedback linearizing 
transformation is 


9 = v k 

Uq = Ai'-X + bisSth 


Eqn 4.27 is inverted to obtain the control-law 


sSth 


- {Vq - A li: x} 


In terms of stability derivatives and pseudo control this is 




1 

2A -Sth 




AqQ AyU AyjUJ ) 


(4.26) 

(4.27) 


(4.28) 


(4.29) 


We now have the freedom to design the signal Ug to provide desired closed loop responses. 
In a model reference setup we design 


"j = 9c + ai(9e-g) + £»o(0 c -0) 


(4.30) 


Analysis of zero— dynamics 

For the longitudinal channel, with a rate command system £ = q. Based on remark 2.1, 
with £ = O=5-q = g = 0 = O, from 4.25 obtain 


A u A 

2A Sth U 2A Sth 


The zero dynamics may be represented by 
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z u 


Au 

7 A \ 

Z M Au 

^ s th 



z w 


7 Af «4 


A *th 


A* 


(4.31) 


(4.32) 
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4.3.1 Analysis of output tracking dynamics 

To analyze the tracking dynamics we start by considering stabilization only, 

9c = Qc = Q c = 0, then 

Vq = -Qi q - a 0 0 (4.33) 

= -kx (4.34) 

Therefore, the symmetric thrust control law becomes 

sS th = £-{-k -A 1>: }x (4.35) 

= — kix (4.36) 

The longitudinal modes in closed loop are then represented by 

A c = A- bk t (4.37) 


4.3.2 Effect of engine dynamics 

Let sy t h be the output of the approximate engine model, 

S±th == ~~K en g SXth -h KengS&thc 

syth = s%th 

and define x' = [ x 3Za ] T Following steps similar to section 4.2.3, 
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x u 

x w 

0 



z tt 

z w 

0 

to 

> 

1 

0 

0 

0 

0 

Ken^ ( n> , - 1 A ^ 
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L 2A Sth + A Q) 
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2 ^ A,u 

2A ^° 

Ix eng 


(4.38) 

(4.39) 


(4.40) 


4.4 Augmentation of Yaw Channel 

As an alternative to providing the pilot with handling qualities in the roll channel, it may 
be easier to provide an autopilot function in terms of heading selection. This is equivalent 
to either a yaw-rate command or heading command, depending on the type of augmentation 
that is desired. The lateral modes are approximated by the linear representation 4.2. where 
a5th represents asymmetric or differential throttle deflection, i.e. the effect of both throttle 
deflections. 

Let A 2i -. represent the second row of A, and 62 the second element of b. A feedback 
linearizing transformation is 


r = Vj. 

Vj. = A 2 , ; x + b 2 a5th 
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(4.41) 

(4.42) 



(4 43) 


Eqn 4.42 is inverted to obtain the control-law 


a5 t h = 


1 


{u f - A 2 , x} 


In terms of stability derivatives and pseudo control this is 

1 


0,5th — 


2Si<„ 


{l/f ~ ftpP ~ ~ 


(4.44) 


We now have the freedom to design the signal &y to provide desired closed loop responses. 
In a model reference setup we design 


= r c + ai(r c — r) + a 0 (^ c - VO 


(4.45) 


4.4.1 Analysis of zero-dynamics 

This analysis is based on linearized expressions 4.2. Neglecting the effects of and Yjth, 
it can be shown that the equations for p and v are in so-called normal form, Eqn2.13. 
Therefore, choose 


V = [pvf 

Prom an analysis similar to section 4.2, with £ = 0 =>• r = 0, obtain 


0,5th — 
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The zero dynamics are represented by 

V = 


r — — 4A o r — T * th o. 
p n*J lp Lv ^ v 

r\ \r *6* , 


Y - V - -^Q 


th 


4.4.2 Closed loop tracking dynamics 

To analyze the zero-dynamics we start by considering stabilization only, 
f c = r c = 7p c = 0, then 

Uf = — r — a 0 ip 

= -kx 

Therefore, the differential thrust control law becomes 

1 


a5 th = 


{-k-4 2 ,,}x 


= — kix 

The lateral modes in the closed loop are then represented by 

A c = A — bkx 


(4.46) 


(4.47) 


(4.48) 


(4.49) 

(4.50) 


(4.51) 

(4.52) 

(4.53) 
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Example 4.4 The closed loop matrix A c for the MD-11 has the following eigenvalues 


-0.780 0.008 ± 0.366* - 0.400 ± 0.301* 


(4.54) 


4.4.3 Effect of Engine Dynamics 


Similar to section 4.2. 

, the effect of the engine dynamics on 

the closed loop stabilization may 

be captured as 
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For the choice of controller gains we may look at a simple root locus analysis where we select 
the ratio K p /K d and vary K d . See the representation in Fig 4.1. The roll channel specifically 
is presented in section 5.4, Fig. 5.3. 
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Table 5.1 Pole Locations in Roll Channel 



open 

loop 

original 

MD-11 

PC , 4 

FBL, k d = 1.0 
k p /k d = 2-0 
ki/k d = 1.0 

dutch roll 

-0.08±0.69i 

-0.291-k0.38i 


spiral & rolling mode 

0.026, -0.80 

-0.13k0.09i 


zero dynamics 


engine 

-2.00 

-2.15 

-2.15 

controller 


-0.35 


FBL poles 



-1.84±0.39i 




-0.55±0.62i 




Figure 5.1: Roll channel root locus, with (left) and without (right) engine dynamics. 


5.3 Pitch rate augmentation 

The zero dynamics are represented by Eqn (4.32), which for the MD-11, consist of two real 
modes with poles at 

-0.1463 and - 0.4581 (5.3) 

The zero dynamics are cancelled, and are therefore ’not visible’, in the feedback linearizing 
control design. The loop dynamics of the linear controller and plant can be approximated 
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as 


(s + kp/kd) 
d ( s* + k d s + k p ) 


(5.4) 


This representation can now be used in a root locus, which will also reveal the effect of the 
engine dynamics. Based on a given ratio kp/kd, the gain kd can be used as the varying gain 
for analysis. Fig 5.2 shows the closed loop results for the nominal MD-11. The left plot 
shows the effect of the feedback linearization, i.e. the open loop modes are cancelled, and in 
its place is a double integrator pole at the origin. The zero dynamics show up as cancelled 
poles. The controller zero at s = —k p /kd draws these integrator poles well into the left half 
plane. The location of the closed loop poles is determined by kd ■ However, the useful choices 
of kd are limited by the engine dynamics. The effect that the engine dynamics have on the 
possible pole locations is indicated in the right plot. Table 5.2 gives the pole locations for 
the original PC A design, and also some examples of pole locations with the current feedback 
linearizing (FBL) setup. 




Figure 5.2: Root loci for the longitudinal MD-11 model, without (left) and with (right) the 
engine dynamics. 
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Table 5.2 Pole Locations in Longitudinal Channel 



open 

loop 

original 

MD-11 

PCA 

FBL, 

k p = 0.36 
k d = 0.72 

FBL, 

k p = 0.26 
k d = 0.83 

FBL, 

kp = 0.50 
k d = 1.20 

phugoid 

-0.01±0.11i 

-0.17±0.19i 

-0.69±0.S2i 

-0.66± 0.58 i 

-0.47±1.09i 

short pd. 

-0.41±0.60i 

-0.44±1.13i 




zero 



-0.45 

-0.45 

-0.45 

dynamics 



-0.15 

-0.15 

-0.15 

engine 

-2.00 

-2.15 

-0.63 

-0.68 

-1.07 


5.4 Yaw channel augmentation 

For the MD-11 [4], the zero dynamics when considering a yaw rate command system axe 
given by Eqn (4.48), which has eigenvalues 


-0.507 - 0.289 (5.5) 

In design of the controller, the effect of the engine dynamics, Eqn 4.55 ,on the tracking 
performance may be captured in a simple root locus analysis as indicated in Fig 5.3. We 
set the ratio a Q /ai = K p jK d and we vary K d . Based on the results we choose a 0 /«i = 
K p /K d — 0.1. In Fig 5.3, K d is varried from 0 — > 0.3. The better response is obtained from 
K d = 0.25 and thus K p — 0.025. The results are poles located at 

-0.184 ± 0.175i, -0.233, -0.223 (5.6) 

This provides a response to a commanded heading with £ = 0.72, u n = 0.25 and T 5 % ze 
16 sec. 
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Figure 5.3: Yaw channel, simple root locus analysis of the effect of engine dynamics. 

6 Stone Soup Simulator 
Implementation 

6.1 Dynamic Inversion for Control Law Design 

In section 3.1 we chose pitch rate/symmetric thrust and roll rate/asymmetric thrust for in- 
put/output linearization with relative degrees of, r = 1 and r = 2 respectively. To implement 
these results using dynamic inversion, we linearize the aircraft dynamics, design controllers 
for the symmetric and asymmetric control signals then transform them to individual engine 
commands. 

6.1.1 Linearized Aircraft Model 

We can expand the rigid body equations of motion using the velocity vector, V = [ u v w] T , 
the angular rate vector, u = {p q r ] T , and the attitude vector, $ = [ $ 6 rp ] T , to make 
up our state vector in the body reference frame [16] 

x = [uvwpqr<f>6 ip] T (6.1) 

The external force and torque vectors in body axis frame are defined 

F = [X Y Z] T (6.2) 
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T = [L M N] t 


(6.3) 


When necessary, the subscript E will be used to identify propulsive (engine) terms. Using 
the following simplifying nomenclature that utilizes some of the McFarland inertia constants 
(q) given in Appendix A. 


T = c 3 L + c 4 N (6.4) 

A = C 7 M (6-5) 

ft = C 4 L + CgN (6.6) 

The force equations are 

u = rv — qw — g' 0 sinO + (X + Xe)/™, (6.7) 

v = pw — ru + g' 0 sin<pcos9 + Y/m (6.8) 

w = qu — pv + g' 0 cos<t>cos9 + Z/m (6.9) 

where g ' 0 is the gravity constant. The moment equations are 


p = c x rq + c 2 pq + T + c 3 L E + c 4 N E (6.10) 

q = c 3 pr - Cg(j? - r 2 ) + A + c 7 M E (6.11) 

f = c 6 pq - c 2 rq + ft + dL E + cgN E (6.12) 

The kinematic equations are 

(f> = p + (qsin(f> + rcos(f>)tan9 (6.13) 

9 = qcos(f> — rsin<t> (6.14) 

ip = {qsincp + rcos<p)/cos9 (6.15) 


The input signals are the symmetric (Me) and asymmetric (N E ) external torques due to 
propulsion. Noting that the symmetric torque component will also contribute a longitudinal 
force component, we can replace Xe with Me/ z t where z t is the engine vertical offset. For 
this class of aircraft, we will further assume a plane of symmetry about the body x,z plane. 
The full nonlinear system of equations are 


’ u ' 


rv — qw — g' 0 sin9 
pw — ru + g' 0 sin<pcos9 
qu — pv + g' 0 cos(pcos9 
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0 
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0. 


\M e ] 

Us. 


(6.16) 


The bar notation indicates normalization by mass of the vehicle (i.e. X = X/m). Using a 
first order Taylor’s series to linearize about a nominal operating point 


x 0 = [no v 0 w 0 po qo r 0 <p 0 9 0 ip 0 ] T 


(6.17) 
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yields the familiar linearized equations of motion 


Aii 
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where the subscripted notation indicates partial derivatives (i.e. ^ = M z ) and the coeffi- 
cients hi are 
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6.1.2 Symmetric Control 

Equations (6.18) and (6.22) show that a symmetric thrust input will result in changes in 
longitudinal angular and translational acceleration. The dominant effect will is a change in 
pitch attitude followed by the secondary effect of a change in longituinal velocities. The focus 
of this section is on the control of pitch attitude and since we are limited to one symmetric 
control effector, it is reserved for pitch control exclusively. Therefore, we chose to ignore 
the effect of symmetric thrust control on aircraft speed. Note, however, that speed stability 
as described in section 4 is inherent to most air-transport aircraft. The linearized pitching 
moment equation, (6.22), confirms that pitch rate/symmetrical thrust is of realtive degree , 
r = 1, since a symmetric thrust change will result in a significant change in the pitch rate. 
Using this equation, the application of Theorem 2.1 is straignt forward. The symmetric 
control transformation is 


A q = Vj (6.24) 

where 


Vg = hi Au 4 - h 2 Av + h 3 Aw + h 4 Ap + h$Aq + h 6 Ar 

+ hjA(j) + /ig A8 + c 7 AM E (6.25) 

The partial feedback linearizing symmetric control law becomes 

A Me = ~-[vq — hi Au — h 2 An — h 3 Aw — h 4 Ap — h$Aq — h 6 Ar 

c 7 

- h 7 A<f> - h 3 A0) (6.26) 

This control law allows use to create the desired response through appropriate construction 
of the psuedo control Vq. 

6.1.3 Asymmetric Control 

Equations (6.21) and (6.23) show that an asymmetric thrust input will result in changes in 
roll rate and yaw rate proportional to the values of c 4 and eg repectively. Because of the 
magnitude of c 4 , trying to implement equation (6.21) would result in a high gain controller 
leading to control saturation issues. For input/output linearization control of roll rate, 
we differentiate (6.21) again. The second derivative of roll rate is significantly effected 
by asymmetrical control through Cg and confirms that roll rate/asymmetrical thrust is of 
approximate relative degree , r = 2. Thus, differentiating equation (6.21) again to obtain the 
second derivative of roll rate 

p = r u Au + T v Av + r w Aw + r p Ap + r,A q + r r Ar ( 6 . 27 ) 

Substituting the results from the linearized equations of motion into the above equation, we 
get the linear approximation for the derivative of the rolling moment equation 

p k kiAu + k 2 Av + k 3 Aw + k 4 Ap + k$Aq + k 6 Ar + k 7 A<j) + k s A6 

+ ( — — + c 7 T q )AM E + (c 4 r„ -f cgT r )ANE (6.28) 
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where the fa coefficients are 
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Again, applying Theorem 2 . 1 , the asymmetric control transformation is 

p = Vp (6.29) 

where 

i/ p = ki Au + £ 2 Aw + A^Au; + A: 4 Ap + k$Aq + A^Ar + k 7 A(f) + faA6 

+ kgAMjz + kigANs (6.30) 

The partial feedback linearizing asymmetric control law becomes 

ANe ~ -r~[i'p — kiAu — fa Aw — fa Aw — A; 4 Ap — A^Ag — faAr 
k io 

- k 7 A(f) - faA 6 - k 9 AM E ] (6.31) 


6.1.4 Engine Control Transformation 

The required change in thrust from the nominal operating point of each engine can be found 
from 
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(6.33) 


A5 2 


1,A Me_ _ AN e , 

2 Z e Ve 


Where z e and y e are the distance of each engine below the nominal C.G. and away from the 
plane of symmetry repectively. Finally, the throttle position for each engine is 

<^1 = <5l(trim) + A5i (6.34) 

S 2 = S2(trim) + A5 2 (6.35) 

The equations developed here can be simplified further by the appropriate choice of the 
equilibrium condition and dropping insignificant terms due to aircraft configuration. 


6.2 Implementation 

Background on the Stone Soup Simulator 

The Stone Soup Simulator (SSS) originated from the Advanced Concepts Flight Simulator 
(ACFS) located at the NASA Ames Crew- Vehicle Systems Research Facility (CVSRF) [17]. 
The ACFS is a full-mission flight simulator representative of a generic commercial transport 
aircraft employing many advanced flight systems as well as features existing in the newest 
aircraft being built today. A miniature version of the ACFS (miniACFS) was developed at 
the NASA Ames Aviation Operations Branch. This desktop flight simulator runs on one 
or more SGI workstations and uses graphical displays to provide ’’soft” interfaces in place 
of hardware controls. The miniACFS provides the flight deck simulation capability for the 
Part-Task Laboratory and the Airspace Operations Laboratory. The Stone Soup Simulator 
(SSS) represents a version of the miniACFS in which proprietary portions of the software 
have been replaced to enable limited distribution. The ’’Stone Soup” name was chosen in 
acknowledgement of the contributions made by the research community which made the 
existence of this tool possible. 


6.2.1 Model Inversion Assumptions 

Because we did not have access to the aircraft characteristics represented by the SSS, we 
used a crude approximation in the aero model for control design. Most notably, the stability 
derivitives used for the dynamic inversion came from the Boeing 747-100 in straight and 
level flight (M=0.2, Sea Level) [18]. This information is included as Appendix B. To keep 
a consistent set of data, the B747 mass properties were also used. Typically, we choose 
the equilibrium condition such that all of the angular velocities, the roll angle and sideward 
component of linear velocity are zero. The coefficients for (6.26) simplify to 
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hi — 
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Then (6.26) becomes 
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MO 

While using the above equations in the SSS, it became apparent that more accurate infor- 
mation is necessary. However, since that information is currently unavailable, we chose to 
adjust the control signals. In particular, since the effectiveness of the propulsive control was 
underestimated, we used scaling to keep their values from saturating when using dynamic 
inversion without the neural network. The resultant scale factors were \ and thf for the 
symmetric and asymmetric controls respectively. The command filters, inversion model and 
neural networks were implemented as described in figure 3.1. For this implementation, we 
assumed that all aircraft states are available for feedback including pitch acceleration (p). 
However, note the latter was not necessary since the pitch acceleration state could have been 
recreated by use of a lead-lag compensator. 

6.2.2 Control Architecture 

Since the SSS aircraft has only the two engines available for propulsion control, only two 
degrees of freedom can be controlled at one time. As stated earlier, speed stability is a 
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inherent in civil aircraft and thus we allow this effect to dictate the speed of the aircraft and 
keep the two control effectors for attitude controls (pitch and roll). The control architecture 
implemented in this study was separated into two parts: the fast inner-loops and the slower 
outer-loops (auto-pilots) as shown in Figure 6.1. 


The inner loop controllers consist of a command filter, a tracking error dynamics com- 



Figure 6.1: Overall Control Archetecture. 

pensator, a linear in the parameters neural network and a dynamic model inversion. The 
command filters are used to prescribe the desired handling qualities to the variables required 
for control, while the tracking error compensator applies these same desired qualities to 
the error signals. The dynamic inversion model then transforms the psuedocontrol to the 
actual control. Since the inversion model is an approximation, the neural network is used 
to reconstruct and cancel model invesion errors. To make resonable comparisons with the 
existing NASA PCA controllers already included in the SSS, simple outer-loop autopilots 
were created for flight path angle and azimuthal (heading) angle. These autopilots have 
much slower time constants than the inner-loops since their primary objective is serve as 
simple pilot models. 

Pitch Channel 

The pitch channel uses a first order command filter followed by a proportional-derivative 
(PD) linear tracking error dynamics compensator. To match the dynamics between these, 
the same natural frequencies and damping ratios are used. Because the command filter is 
only first order, r = 2(u; n was used to approximate the relationship of the real (first order) 
and complex conjugate (second order) poles. The linear in the parameters neural network 
was implemented as described in Reference [3] to recontrust the model inversion errors. The 
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neural networks employ the following scaled input signals: 


V = (V — 300)/500 
H = (H - 30000)/40000 
i - e~ Vp 
Up ~ 1 + e - ** 

1 - e~ v « 


(6.38) 


where V = y/u 2 + 1> 2 + 1 c 2 and H is the altitude. In addition to the above signals, the NNs 
will require some of the states {(f>,9,p,q) and some of the command filter signals (g c ,p c ,p c ). 
The input vectors, described in Reference [3], for the pitch channel are constructed from a 
0.1 bias and the above signals as follows 

Cb = [0.1 V V 2 H] t 

c 2 = [0.1 Pg 9 <t> q q c ] T 

C 3 = [0.1 a p] T (6.39) 


The outer-loop for the pitch channel is a simple flight path angle hold autopilot consisting of a 
simple integrator. Since we are exciting the phugoid mode to control pitch, the assumption 
that a is negligible is warrented. The flying qualities of this controller is controlled by 
specifying r g . The command filter, tracking error compensator and flight path angle hold 
auto pilot are shown in Figure 6.2. Control saturation is avoided by limiting the error signal, 
e 7 = 7 cmd ~ 7, by the equivalent error produced with a maximum pitch rate, 


^7 — QlimitTg 


(6.40) 


Roll Channel 

The roll channel differs from the pitch channel in a number of ways. The primary deviation 
is due to the relative degree of the control in this channel. As stated earlier, the relative 
degree of the roll control is two and thus the psuedo control is related to p. In order to 
provide this information, the command filter needs to be at least second order. However, to 
soften acceleration spikes, a third order command filter with a fast pole at s = —aC^n is 
used. The associated transfer function is 


Pc 

Pcmd(s) (s + a(uj n )(s 2 + 2 (ui n s + a;^) 


(6.41) 


The command filter is followed by a proportional-integral (PI) linear tracking error dynamics 
compensator. Again, to match the dynamics between these, the same natural frequencies 
and damping ratios are used in the PI compensator and the dominant poles of the command 
filter. Using the scaled inputs descibed for the NN in the pitch channel and the appropriate 
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Figure 6.2: Pitch channel a) Command filter, b) Tracking error compensator, c) Flight path 
angle hold autopilot. 
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sigmoidal function, the input vectors for the roll channel NN are constructed as follows 

Cj = [0.1 V V 2 H) t 

C 2 = [0.1 u p <p q c q c ] T 

C 3 = [0.1 a P] T (6.42) 

Note that Ci and C 3 for the pitch and roll channels are identical. The outer-loop for the 
roll channel is an azimuthal angle (heading) hold autopilot consisting of a cascade of two 
simple integrators on ip and <p. The gain ^ is required to model this channel as two 
integrators, where V# = Vcos8 is the horizontal plane aircraft velocity. Figure 6.3 shows 
the block diagrams of the roll channel command filter, tracking error compensator and the 
heading hold autopilot. The flying qualities of beading hold autopilot is straight forward. 
First, consider the simplified representation of this loop with the aircraft replaced by two 
integrators as shown in Figure 6.4. The transfer function is 


£(g) = KpK* 

^cmd(s) s 2 + Kqs + K^K,p 


(6.43) 


with 2(u n = K# and w 2 = K^K^. This allows us to specify the flying qualities in terms 
of the desired natural frequency and damping ratios. Again, this loop, by design, is made 
slower than the inner loop. Control saturation was avoided by limiting the error signal, 
= iPcmd — ip, by an equivalent error for a maximum roll angle or a maximum roll rate 
propagated through the appropriate gains, 

e* = rnin(<p limit , + <p) g (6.44) 


6.2.3 Simulation Results 

The SSS was initiated in the nominal aircraft configuration and at the design operating con- 
dition (Vr=180kts, ff o6 ,=2000ft, 5f laps =0°, slats retracted, CG=0.2875ca//ic) and allowed to 
reach a steady state. The conventional flight control surfaces were failed and, immediately 
following, the Intelligent Flight Control System (IFCS) was engaged. A pulsed angular rate 
doublet was then sent to the command filters which generated the desired command signals 
based on the predetermined flying qulaities. Note that separate runs were performed for the 
roll rate and pitch rate doublets to show the rate command attitude hold responses of the 
inner-loops without the outer-loop controllers. 

Figure 6.5 shows the aircraft response using the model inversion only and also with neu- 
ral network adaptive augmentation for a pitch rate doublet command. The command filter 
frequency and damping ratio were chosen to be 0.5 rads/sec and 0.8 respectively. The mag- 
nitude of the pitch rate command (about 0.25 deg/s) was selected to represent a nominal 
transport aircraft maneuvering rate while the pulse width keeps the aircraft from depart- 
ing too far from the initial conditions. Figure 6.6 shows similar responses for the roll rate 
doublet command. Here, the command filter frequency and damping ratio were chosen to 
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Figure 6.3: Roll channel a) Command filter, b) Tracking error compensator, c) Flight path 
angle hold autopilot. 
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Figure 6.4: Simplified roll channel autopilot loop. 

be 1.5 rads/sec and 0.8 respectively. The magnitude of the roll rate command is about 0.60 
deg/s and is far below the nominal rate for a transport aircraft. This, in part, is due to 
the reduced roll control authority of the asymmetric engine inputs. Again, the pulse width 
was selected to keep the aircraft in a ’’proper” attitude. These plots show that the neural 
network is compensating for the dynamic inversion error due to the crude approximation 
used. 

To show the stability and versatility of this control architecture, runs were made with an 
agressive simultaneous pitch rate and roll rate doublet for the design point as wall as ’’off 
design” points. The first case was to alter the aircraft configuration by lowering the flaps 
to Sp = 20° and extending the slates to 5s = 20° at the design flight condition. The second 
case is the nominal aircraft configuration at the cruise condition (Vr=285kts, /f O (, s =30000ft). 
Figure 6.7 shows the baseline aircraft response to this manuever at the nominal design point. 
Figure 6.8 shows the response for the aircraft configured with the flaps and slats extended at 
the same flight condition. Figure 6.9 shows the baseline aircraft response for a cruise condi- 
tion (30,000ft altitude, Mach 0.7). Note that the magnitudes of the doublets were halved to 
keep the aircraft in a ’’proper” attitude during this high speed manuever. The oscillations 
exibited in the roll rate response is due to unmodeled nonlinearities in the engine dynamics, 
in particular, rate saturation. Figure 6.10 is the response to a changing aircraft configuration 
(i.e. 5p = 0°to20°) while the inner-loops regulate zero pitch and roll rates compared to the 
dynamic model inversion only case. All these figures indicate that the neural network will 
compensate for aircraft configuration and flight condition changes across the flight envelope 
relative to the design point for the model. 

To compare responses to the SSS NASA PCA, the outer-loops were engaged to hold aircraft 
flight path angle and heading during a flap deployment. The result is shown in Figure 6.11. 
The saturated control lead to integrator wind up that caused the large overshoot in pitch rate 
at approximately 30 seconds. Figures 6.12 and 6.13 show the vehicle responses to step flight 
path angle and heading commands respectively for the nominal design point also compared 
to the NASA PCA. Lastly, Figures 6.14 and 6.15 show the same for the up and away (cruise) 
condition. These figures show that the neural network architecture and the associated outer- 
loops provided very similar handling qualities in flight path angle and heading responses as 
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the NASA’s SSS PCA for both fight regimes. However, the pitch rate and roll rate responses 
were quite different. The NASA PCA displayed more oscillations in the pitch rate and larger 
variances in roll rate to achieve similar results in flight path angle and heading respectively. 
It should be noted that the rate limits (like the inner-loop responses) and bank angle limit 
were reduced to accommodate the high speed case. 
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Figure 6.5: Inner-Loop Rate Command Attitude Hold Response in Pitch Channel 

(V=180kts, H=2, 000ft, 5 F =0deg). 
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Figure 6.6: Inner-Loop Rate Command Attitude Hold Response in Roll Channel (V=180kts, 
H=2, 000ft, <S F =0deg). 
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Figure 6.7: Simultaneous Roll Rate and Pitch Rate Doublets (V=180kts, H=2, 000ft, 
5 F =0deg). 
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Figure 6.8: Simultaneous Roll Rate and Pitch Rate Doublets (V=180kts, H=2, 000ft, 
5F=20deg). 
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Figure 6.9: Simultaneous Roll Rate and Pitch Rate Doublets (V=285kts, H=30, 000ft, 

<*F=0deg). 
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Figure 6.10: Flap Deployment with Pitch Rate Regulation (V=180kts, H=2, 000ft, <S F =0 to 
20deg). 
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Figure 6.11: Flap Deployment with Flight Path Angle Regulation (V=180kts, H=2, 000ft, 
Sp = 0 to 20deg). 
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Figure 6.12: Outer-Loop Flight Path Angle Step Response (V=180kts, H=2, 000ft, 5F=0deg). 
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Figure 6.14: Outer-Loop Flight Path Angle Step Response (V=285kts, H=30, 000ft, 

<J F =0deg). 
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Figure 6.15: Outer-Loop Azimuthal (Heading) Angle Step Response (V=285kts, H=30,000ft, 
<SF=0deg). 
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7 Recommended Further 
Developments 

This section provides a brief overview of what we envision as the major aspects of further 
investigation. 

Model accuracy and data availability 

The implementation of the controller architecture was hampered by limited availability of 
data on the SSS aircraft model, as well as by limited access to the known data. The crude 
knowledge of the model still allowed for a good demonstration of the flexibility of the con- 
troller architecture, demonstrating the fact that it does not rely on detailed aerodynamic 
and performance data as required in more classic gain scheduled approaches. However, pru- 
dence in the design of flight controllers, as well as the desire for clean implementations, 
suggests that more accurate knowledge of the SSS aircraft model should be used in further 
development. 

Single-Hidden-Layer NN implementation 

The current demonstration was based on the controller architecture augmented with a linear 
in the parameters NN. Linear and non-linear single-hidden-layer (SHL) NNs fit within the 
same controller framework. The SHL-NNs are more powerful in their approximation prop- 
erties, and are therefore preferred for highly nonlinear mapping tasks. In nominal operation, 
it is hard to demonstrate their advantage over linear NNs. However, the cases we wish 
to consider include the possibility of asymmetric failures, and a variety of possible control 
effectors. Therefore, the nonlinearities may become significant, and the use of SHL-NN is 
preferable. 

Tolerance to unanticipated failures 

As more detailed knowledge of the SSS aircraft model becomes available, we wish to evolve 
the controller into a fully integrated IFPCS. This controller will enhance the flight safety 
in case of unanticipated critical failures by utilizing a variety of conventional aerodynamic 
surfaces and configuration controls. For example, hydraulic surfaces, speed brakes, electro- 
mechanically driven spoilers, flaps and or slats deployment. The controller architecture as 
presented, with access to a variety of these effectors, can provide tolerance to unanticipated 
failures. The implementation of this work will require access to the various dynamic models 
of the configuration surfaces within the SSS code. 

Taking advantage of the maximum available propulsion performance 

The results in section 5 include oscillatory behaviour which may be caused by interaction of 
engine dynamics and the inverting controller. This demonstrates an inherent drawback of 
feedback linearizing controllers. An illustration of this effect is given in Fig 7. The oscillatory 
interference would usually be avoided by further bandwidth reduction of the controller. In the 
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Potential of Direct Adaptive NN Output Feedback Control 


i, a 



The effects of luunodeled actuator 
dynamics on inverting controllers 



Dealing with the effects of actuator 
dynamics witliout altering its time 
constants, yet taking advantage of its 
maximum performance. 



Experience in application to unmodeled 
control rotor, is representative of unmodeled 
engine dynamics in IFPCS application 


Figure 7.1: Illustration of the effect of unmodeled actuator dynamics. The effect 
of the unmodelled control rotor in the high bandwidth control of the R-50 
autonomous helicopter, is similar to the effect of unmodeled engine dynamics 
for the IFPCS. 
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IFPCS this is undesirable because of the performance requirements for safety. Furthermore, 
one would like to take advantage of the maximum available propulsion performance, without 
attempting to alter its inherent dynamics. 

Preliminary results, based on a recent development in Direct Output Feedback Control 
Ref [19, 20, 21], demonstrate the ability to provide robustness with respect to unmodeled 
actuator dynamics and delays. The IFPCS is an ideal platform for application of this recent 
development. 

Dealing with actuator saturation 

Stability and robustness with respect to actuator saturations are subjects of many recent 
investigations, e.g. Ref [22], A controller structure similar to that of section 3 is possible, 
using a form of predictive control Ref [23], and allowing for both rate and position saturations 
Ref [24]. Rather than a pure inverting controller, this would provide for an optimal controller. 
The optimal controller performance is guaranteed in the face of actuator saturations. When 
not saturated, the controller signal is very similar to that of the inverting controller, with 
the additional effect of control weighting. The optimal controller may be augmented with an 
adaptive NN in a similar fashion as presented in this report. This would lead to an adaptive 
IFPCS with some robustness with respect to actuator saturations. 
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8 Conclusion 

This document reports on the first phase of the development of an Intelligent Flight Propul- 
sion Control System (IFPCS). The approximate Feedback Linearizing (FBL) controller of 
previous work of the authors was integrated with the approximate Input/Output Lineariza- 
tion formalism from recent literature. A partial FBL controller, based on a crude approxima- 
tion of the actual aircraft, is developed. Adaptive Neural Networks (ANNs) are introduced 
to compensate for the inversion error caused by cascading the linear inversion model with 
the actual aircraft. The resulting controller design is used to drive the throttle inputs of 
both wing mounted engines to achieve desired pitch rate and roll rate tracking. The actual 
aircraft is represented by the NASA Stone Soup Simulator representing a generic midsized 
transport aircraft. 

Results indicate that, despite using only crude knowledge of the aircraft, performance 
similar to extensively gain scheduled design is achievable. The addition of ANNs provides 
robustness with respect to the uncertainties and unmodeled nonlinearities in the aircraft. 
Therefore, the nonlinear adaptive design enhances the conventional controller in significant 
ways. Specifically with respect to, 

• consistency of response over a variety of failures, 

• applicability to different aircraft without redesign, 

• reduction in extensive gain scheduling research & design. 

Immediate further work will include a demonstration of propulsion only control in some 
unanticipated (including asymmetric) failures. 

The architecture, as presented in this report, provides a powerful bases for an IFPCS as 
envisioned by the NTSB & NASA administrators. Specifically; 

• the approximate I/O linearizing setup can be extended to output feedback designs, 

• the application can be enhanced by utilization of all available primary and secondary 
control surfaces, with straight forward adjustments for different bandwidth capabilities 
and control authority limits. 

A continuation of the research & development will include: 

• Improvement of the inversion model parameters when more data becomes available. 

• Implementation and demonstration with the single-hidden-layer neural network. 

• Development of a fully integrated IFPCS, which will make use of spoilers, flaps and 
slats, to investigate to what extend handling qualities may be maintained under various 
unanticipated failure scenarios. 

• Use of recent developments in direct output feedback, which is not limited by the need for 
bandwidth separation with engine dynamics, and thus can take maximum advantage 
of available propulsion. 

• The augmentation of the existing controller architecture to allow for actuator and 
engine rate and position saturations. 
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McFarland Inertia Coefficients 
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B Data for the Boeing BT47-400 



X(lb) 

Y(lb) 

Z {lb) 

u{ft/s) 

-366.1 

0 

-3538 

v{ft/s) 

0 

-1559 

0 

w{ft/s) 

2137 

0 

-8969 

w(ft/s 2 ) 

0 

0 

585.1 

p(rad/s) 

0 

0 

0 

q(rad/s) 

0 

0 

— 1.090e5 

r(rad/s) 

0 

0 

0 


Weight = 



L(ft * lb) 

M(ft * lb) 

N(ft * lb) 

0 

3779 

0 

-86120 

0 

39750 

0 

-57170 

0 

0 

-7946 

0 

-1.370 * 10 7 

0 

— 6.688e6 

0 

— 1.153e7 

0 

4.832e6 

0 

— 1.014e7 

5.640 * 10 b slug * ft 2 


1.42 * 10 7 slug 

* ft 2 


3.23 * 10 7 slug 

* ft 2 


4.54 * 10 7 slug 

* ft 2 


8.70 * 10 5 slug 

*ft 2 
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